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ABSTRACT 

We report an analysis of the dynamical structure of clusters of galaxies from a survey of photometric 
and spectroscopic observations in the fields of southern Abell Clusters. We analyze the galaxy velocity 
field in extended regions up to lhr x Mpc from cluster centers and we estimate mean velocity dispersions 
and their radial dependence. Only one from a total number of 41 Abell clusters does not correspond to 
a dynamically bound system. However, four of these bound objects are double clusters. We estimate 
that 20 % (7 clusters) of the 35 remaining are subject to serious projection effects. Normalizing the 
clustercentric distances by means of the overdensity radius r2oo, and the velocity dispersion profiles 
(VDPs) by the corresponding mean cluster velocity dispersion, we computed the average VDP. Our 
results indicate a flat behavior of the mean VDP at large distances from the cluster center. Nevertheless, 
we found that for the inner part of the clusters (r/r2oo < 1) the VDP is up to a 10% smaller than at 
larger radii. 

Subject headings: galaxies: clusters: general, individual — surveys 



galaxy systems grow by aggregation of smaller structures 
formed earlier. Therefore, we expect a significant degree of 
substructure in clusters of galaxies if the remnants of the 
accretion of groups in the recent past has not been erased 
by the dynamical relaxation of the clusters. The substruc- 
ture in rich clusters has been extensively analyzed in re- 
cent years (Dressier & Shectman (1988), West & Bothun 
(1990), Zabludoff et al. (1993), Girardi et al. (1997), 
Solanes et al. (1999)). The results are consistent with sub- 
structure in most of the cases studied, irrespective of the 
samples and method of analyses used. West & Bothun, 
(1990) made an analysis of substructure in clusters of 
galaxies and their surroundings. The authors developed a 
technique that is sensitive to correlations between galaxy 
positions and local kinematics finding little evidence for 
substructure in the inner regions, and significant depar- 
tures from a relaxed substructure-free systems in the ex- 
ternal regions. More recently, Biviano et al. 2002 realized 
a detailed analysis of the consequences of substructure on 
luminosity and morphological segregation. These authors 
find that the number of galaxies in substructures decrease 
markedly toward the cluster center and report differences 
in the properties of galaxies depending whether they be- 
long to substructures or not. These differences are also 
present in the the dynamical properties of galaxies. 

Escalera et al. (1994) provide an extensive discussion 
of the presence of substructure in clusters of galaxies by 
using galaxy positions and redshifts. In their studies a 
multi-scale analysis is adopted that considers the kinemat- 
ics as well as the wavelet transform providing estimators 
of the degree of substructure. Other works (see for in- 
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1. INTRODUCTION 

Analysis of large scale structure formation may greatly 
benefit from studies of the dynamics of clusters of galax- 
ies. Measurements of galaxy velocity dispersions in clus- 
ters provide reliable estimates of cluster masses and a di- 
rect normalization of the primordial mass power spectrum 
(Eke, Cole & Frenk, 1996). Moreover, the velocity field in 
the extended halos of clusters may set additional impor- 
tant constraints to the formation of structure as well on 
the mean density parameter of the universe. 

There have been several recent studies on the dynam- 
ics of clusters of galaxies, see for instance Girardi et 
al (1993), Zabludoff et al (1993), Collins et al (1995), 
Mazure et al (1996), Fadda et al. (1996) Alonso et al. 
(1999). The resulting distribution function of velocity 
dispersions from the ESO Nearby Abell Cluster Survey 
(ENACS) given by Mazure et al (1996) is in agreement 
with the distribution of cluster X-ray temperatures, sug- 
gesting /3 = cr/iTOft/ '(kTx) — 1. The velocity dispersion 
profiles (hereafter VDP) may provide a useful tool for the 
study of the dynamics of clusters of galaxies. The anal- 
ysis by Fadda et al. 1996 is consistent with a tendency 
of flat VDP in rich Abell clusters. Jing and Borner 1996 
investigated the VDPs of clusters for several cosmological 
models. They found that on average, VDPs decrease with 
the cluster radius in every model up to lh~ l Mpc from the 
cluster center. Also, these authors found that the slope of 
the profiles are different in different models, being steeper 
in lower-f2 models than in higher-fi models. 

In the hierarchical scenario of structure formation, 
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stance Fadda et al. 1996) consider velocity gradients and 
anisotropy of galaxy orbits. Extensions of the different 
methods of analysis can provide new useful quantitative es- 
timates of substructure, essential for a better understand- 
ing of the dynamics of clusters of galaxies. 

The dynamics of clusters of galaxies may also be stud- 
ied from information in the X-ray band. X-ray emission 
detected in a large fraction of clusters of galaxies pro- 
vides an invaluable observational material. Several prop- 
erties of the clusters and the intra-cluster medium may 
be addressed with this information. For example, the 
global mass distribution, the dynamical state and the evo- 
lution with redshift , the composition of the intra-cluster 
medium, etc. White (1999) presents an elegant method- 
ology to recover the spatial properties of the intra-cluster 
gas from X-ray observations. From the deconvolution of 
ASCA satellite X-ray data, he finds a large fraction (90%) 
of clusters consistent with isothermality. These results 
are in conflict with the Markevitch et al. (1998) analy- 
sis from a sample of 30 clusters where most show steeply- 
declining intra-cluster temperature profiles. In their anal- 
ysis of ASCA resolved spectroscopic data these authors 
obtained projected temperature profiles and in many cases 
two-dimensional temperature maps concluding that the 
gas temperature varies within a factor 1.3-2 or greater 
within the clusters. 

The conflicting evidence for isothermality of the intra- 
cluster medium show that the information on the VDP 
for clusters may add important information to the sub- 
ject. On the other hand, the tendency for subclustering 
to occur at large distances from cluster centers encour- 
ages us to explore the outer regions of clusters of galaxies. 
In this paper, we analyze the radial velocity distribution 
in regions extending up to 7/i~ 1 Mpc in projection from 
the Abell cluster center (we adopt Ho—lOOhkms^ 1 and 
qo =0.5). We provide a detailed analysis of each individ- 
ual cluster providing the degree of substructure and an 
estimate of the VDP at large distances from the cluster 
center. In section 2 we describe the method of analysis 
of substructure. Section 3 deals with the identification of 
the clusters and the projection effects which significantly 
affect the measurements of velocity dispersions. Section 
5 provides the estimates of mean velocity dispersion and 
the correlation with richness counts as well as the velocity 
dispersion profile of several clusters. 

2. DATA 

The SARS survey (Southern Abell clusters Redshift Sur- 
vey, Way et al. 2002) comprises Abell 1958 and Abell, 
Corwin & Ollowin 1989 (hereafter ACO) clusters R > 1, 
principally in the region < 5 >-65 and 5h> a < 21h 
(avoiding the LMC and SMC), with b < -40. Galaxies 
were selected from the APM catalog (Maddox el al. 1990). 
Galaxies brighter than jur = 19 and with surface bright- 
ness within 1.5 x 1.5 deg 2 centered on the cluster were 
pre-selected. Target galaxies were selected at random and 
the final completeness is roughly constant up to an appar- 
ent magnitude ~ 18 and it is of the order of 75%. 

The observations were carried out with the 2.5m 
DuPont telescope at the Las Campanas Observatory, 
Chile. The multi-fiber spectrograph (Shectman 1989) was 
used. Fibers are connected to a Boiler & Chivens spec- 
trograph attached to a 2DFrutti detector (2DF). The un- 



known spectra were calibrated using software packages 
within IRAF, following essentially the method described 
by Way et al. (1997). 

From the wavelength calibrated spectra was obtained 
the respective radial velocities cz of the unknown spec- 
tra by using two different, independent methods: i) The 
Fourier Cross Correlation Technique, where two fouricr 
transformed spectra, the unknown object and a known 
template, are multiplied together to obtain the fourier 
transform of their correlation function (with RVSAO, 
Tonry et al. 1979). ii) "By eye" identification of absorp- 
tion lines and computation of cz with the task rvidlines. 
The final sample consists in more than 4000 galaxies with 
redshift estimates in 41 clusters. Mean cluster's redshift 
run from 0.06 to 0.16 with a mean around 0.088. 

3. ALGORITHM FOR SUBSTRUCTURE DETECTION 

We have applied two different techniques in order to de- 
tect substructure in clusters. These techniques are com- 
plementary in the sense that they are mainly designed to 
remove large structures along the line of sight, and smaller 
systems in three dimensions. Many clusters present dou- 
ble structures in the redshift distribution (e.g late stage 
of a cluster-cluster merger). Ashman et al. (1994) dis- 
cuss a statistical technique for detecting and quantifying 
bimodality known as mixture modeling or the KMM al- 
gorithm. The scheme is based on the application of algo- 
rithms that fit a certain number of substructures in red- 
shift space to the data and it is determined the best fitting 
model. This technique is the base of commonly adopted 
procedures used to analyze astronomical data sets and it 
assesses the statistical significance of bimodality providing 
objective ways of dividing the data into sub-populations. 
As discussed by Ashman, Bird & Zepf (1994) the KMM 
technique has broad applicability in the analysis of as- 
tronomical data. We have applied this technique to the 
redshift distribution regardless the angular position of the 
galaxies in the field of the cluster. Based on a preliminary 
inspection of the data we propose the number of struc- 
tures with their corresponding mean radial velocities and 
velocity dispersion that approximately represent the red- 
shift distribution around the cluster. This procedure is 
restricted only to those structures with overlapping red- 
shift distributions. Then we apply the KMM algorithm 
and consider a multiple peak structure when the confi- 
dence level to have proposed model is bigger than 90% 
(for details see Ashman et al. 1994). We have consid- 
ered different possibilities: i) when the proposed model 
for multiple-peak-structure is rejected we consider that the 
redshift distribution corresponds to a single cluster; ii) if 
the confidence level of the proposed model is bigger than 
the 90% and at least 70% of the galaxies belong to the 
same structure we assume a single cluster and discard the 
outlying groups that will deserve a detailed study in a fu- 
ture work.; iii) when most of the galaxies belong to two 
separate structures of similar sizes we assume the presence 
of two clusters and we perform the corresponding analy- 
sis. It should be noted that the above technique works 
properly when the sub-structures are representative of an 
important number of galaxies. In this work only structures 
with at least 10 galaxies are considered. 

In the hierarchical model for large scale structure for- 
mation, clusters of galaxies are the result of a continu- 
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ous process of accretion of small structures like groups of 
galaxies. Therefore, a considerable number of galaxies are 
expected to be found around clusters that arc not bound 
to the main system and therefore will bias the velocity dis- 
persion estimate if they are included in the analysis. This 
problem is particularly serious if large distances from the 
cluster center are considered as is the case in the present 
work. Some of these groups of galaxies can be located at 
a similar redshift to the cluster, therefore, they are very 
difficult to detect in the redshift space . 

In our data set, and for each cluster, we analyze the real 
nature of visually identified group candidates with a tech- 
nique that considers both the projected position and red- 
shift of the galaxies using the compactness of the projected 
distribution and the departure from the mean dynamical 
properties of the cluster. Three different parameters are 
used: 

i) a S parameter similar to that defined by Dressier and 
Shectman 1988: 



ber of galaxies with measured redshift in the line of sight 
of the cluster (N to t)', the number of galaxies assigned to 
the cluster (N c i u ) and the cluster mean radial velocities. 



S 2 = (U/a 2 )[(v local - vf + (a local - a) 2 } 

where a and ai oca i are the velocity dispersion of the clus- 
ter and the proposed group structure respectively and v 
and vi oca i are the corresponding mean velocities. 

ii) A parameter C, that provides a measure of the com- 
pactness of groups and is computed as: C =< dnni oca i > 
I < dnn group > where < dnn group > is the average pro- 
jected distance of the nearest neighbors to each of the ng 
members of the proposed group and < dnni oca i > is com- 
puted in the same way but for the nearest ng galaxies in 
the neighborhood of the proposed group. 

iii) An isolation parameter / = dnng/ < dnn group > 
where dnng is the distance to the nearest neighbor galaxy 
to the group. 

We compute the variable G = S + C + I and we calculate 
the mean (G) and the dispersion gq for each cluster. A 
given group candidate is to be removed if the value of G 
for the group is at least 2 standard deviations, 2ctg, away 
from the cluster mean value (G). The adopted threshold 
is the result of Monte Carlo simulations which show that 
this threshold is suitable to remove structures. For the 
five most regular clusters in our sample we reassigned the 
polar angle of every galaxy with respect to the cluster cen- 
ter. This procedure removes group structures and leaves 
unchanged the radial galaxy density profile of the clus- 
ter. For the mock data, we compute G identifying mock 
groups of galaxies finding that none of these chance groups 
have G-(G) > 2a q (with (G) and oq computed from the 
original data). 

In spite of the fact that our sample of galaxies in clus- 
ters is not magnitude limited the above procedure should 
give no biased results provided the galaxies are randomly 
selected from a complete sample. 

4. ANALYSIS 

As a result of the two techniques described above we 
have removed 19 structures in 14 clusters from our total 
sample of 41 Abell clusters. We find that 4 Abell clus- 
ters appear as two different systems in redshift space, and 
one Abell cluster is completely spurious. In the following 
sections we discuss different properties of the resulting 44 
clusters. Table 1 shows the Abell number; the total num- 



4.1. Cluster identification 

In our analysis we have only considered clusters selected 
by Abell 1958 and ACO. Several authors (van Haarlem et 
al. 1997 and references therein) have discussed the conse- 
quences of the projection effects when clusters are selected 
in 2-D catalog. Redshift surveys provide precise informa- 
tion on the reality of the clusters selected. As a result 
of our analysis we found only one spurious cluster (Abell 
3159) while the rest appear as real concentrations in red- 
shift space. Nevertheless, 11 clusters present more than 
one concentration in the redshift space, thus in projection 
they appear as richer clusters. Of the total of 40 Abell 
clusters 28 appear as a single concentration in the total 
redshift range, while the rest have been systematically en- 
hanced by projection effects. We consider that a cluster 
is significantly affected by projections when the number 
of galaxies in groups or other cluster-like structures along 
the line of sight is comparable to the number of confirmed 
cluster members. Besides the projection effects, and af- 
ter the removal of structures previously described, several 
clusters present significant evidence for substructure on 
different levels. This substructure can affect the analysis 
of the dynamics of clusters. In particular, the estimate 
of the velocity dispersion may be significantly affected by 
substructure. 

We will use the term "relaxed" cluster to describe 
a system that is free from substructure with a single 
nearly gaussian redshift distribution after the subtraction 
of structures using the procedure discussed in Section 3. 
This classification will be used to define sub-samples of 
clusters. Our original cluster sample is not volume com- 
plete, and the above definition is used to select subsets of 
clusters in order to cross-correlate general properties such 
as richness counts, mean velocity dispersions, etc. Figures 
la and lb show the mean redshift distribution of the total 
sample clusters and those classified as "relaxed" respec- 
tively where it can be appreciated the similarity of the 
redshift distributions. 
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Figure 1. Histogram of mean radial velocities of the clusters: 
a) the total sample; b) "relaxed" clusters. 

4.2. Substructure properties 

We have analyzed the properties of the different struc- 
tures removed from clusters. As a result of the algo- 
rithm of group detection we find an average velocity dis- 
persion <jg = 295 ± 180km/ sec and a mean extension 
D = 0.44 ± 0.28/i _1 Mpc. These structures which com- 
prise 12% of the total number of galaxies in the clusters 
have, on average, a difference of mean velocity respect to 
the parent cluster AV — 921 ± 393km/sec. Our values 
of a and D are typical of group of galaxies. Neverthe- 
less, the average extension of our groups are larger than 
those derived by Girardi et al. 1997 (~ 0.2h~ 1 Mpc). The 
KMM technique for substructure detection tends to iden- 
tify systems at larger distances from the cluster center 
(AV = 1515 ± 304km/ sec). Nevertheless, the mean ve- 
locity dispersion of these structures (288 ± 128km/ sec) is 
similar to the a g derived by the group detection algorithm. 

4.3. Individual objects 

Several clusters in our analysis deserve individual atten- 
tion due to different peculiarities of their properties. 

The rcdshift distribution along the line of sight of A2819, 
A2871, A3107 and A3223 show two similar structures not 
physically connected. In Table 1 each of these clusters are 
named with the original Abell number plus an "a" and "b" 
respectively. 

As an example of the application of the method we com- 
ment on the cluster A2734 which presents a double peaked 
structure in redshift space. The smallest peak has approx- 
imately half the number of members of the main structure. 
It was removed since the probability to have two different 
structures is 99 %. 

A380 presents some evidence of a double structure with 
a mean velocity difference of 1407 km/s. The probability 
to have two different structures is larger than 90%. Nev- 



ertheless, the low number of galaxies (25) involved in our 
analysis introduce some doubts into our conclusions. The 
values quoted in Table 1 correspond to a single cluster. 
Assuming two different structures we find the following 
values: (V) = 31440 km/s, cr=408 for the nearest struc- 
ture (14 galaxies) and < V >= 32847 km/s, cr=314 for 
the second (11 members). 

Besides the clusters that appear as double in the red- 
shift space, A380, A487, A2915, A3142, A3153, A3844, 
A3864 present strong projection effects due to the pres- 
ence of several structures such as groups of galaxies along 
the line of sight. 

A3111: This cluster shows some evidence of large scale 
substructures in the plane of the sky. Our algorithm does 
not work properly for this type of substructure, there- 
fore, the cluster was taken as a single structure and the 
value quoted in Table 1 (943 km/s) could be biased high. 
Nevertheless, the discrepancy with the velocity dispersion 
derived by Fadda et al. 1996 (159 km/s) can not be ex- 
plained. If we arbitrarily restrict ourselves to the cen- 
tral region of the cluster (up to 2.5 Mpc/h in diameter) 
where no evidence of substructure is present we derive a 
a = 734fcm/s. This value must be taken as a lower limit 
for the mean velocity dispersion of A3111. The value de- 
rived by Fadda et al. 1996 is probably biased by the low 
number of confirmed members in their sample (12 galax- 
ies) while our analysis is based on more than 50 cluster 
members. 

A3151 presents a group of galaxies in the very center of 
the cluster with a mean velocity differing by more than 900 
km/s with respect to the main cluster. This is nearly the 
same difference between our estimate of the cluster mean 
velocity and the value derived by Fadda et al. 1996. Their 
estimate is the result of 14 galaxies, and by chance they 
selected galaxies from the projected group instead of the 
main cluster. 

A3223 appear as two separate structures in the plane 
of the sky and hence was treated as two different clus- 
ters. These two clusters also show important differences 
in their dynamical properties. The second concentration 
was identified by the APM selection criteria and is named 
as APMCC 479. 

After the removal of groups and besides those indicating 
double structures A1750, A3111 A3135, A3764 A3915 still 
present some evidence for substructure in redshift space 
or in the plane of the sky. 

A3159: The redshift distribution in the line of sight of 
this cluster shows the presence of several groups, never- 
theless, none of these groups can by classified as a cluster. 
We suggest this system is a spurious cluster identification. 

A2778 and A3153 are two clusters poorly defined both 
in the plane of the sky and in the redshift space where 
the presence of gaps suggest the possibility of substruc- 
ture. More redshifts are needed in order clearly under- 
stand these clusters. The values of a quoted in Table 1 for 
these two clusters must be taken with caution, especially 
in the case of A3153 where the redshift distribution can 
be also consistent with several groups instead of a single 
cluster. 

5. VELOCITY DISPERSION ESTIMATES 

After the redefinition of structures as defined in Section 
2 we have computed the mean velocity dispersion for each 
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cluster. Based on the ROSTAT routine (see Beers et al. 
1990) we have used robust mean and scale estimators. We 
have applied relativistic corrections and have taken into 
account velocity errors. Considering the typical number 
of redshift confirmed cluster members (usually > 20) we 
have considered the biweight estimate for both the mean 
cluster radial velocity and the velocity dispersion. Errors 
are based on the statistical jacknife. The derived values 
are shown in Table 1. Figure 2a shows the values of the 
mean velocity dispersion for clusters in the range 200 - 
1100 km/s with a mean w 600 km/s indicating that in 
our study we have included low mass systems (probably 
groups) as well as massive clusters of galaxies. Figure 2b 
shows the same distribution but only for those clusters 
classified as "relaxed". As it can be appreciated no dif- 
ferences are present between both sets of data indicating 
that contamination by projection effects are seen at some 
degree in all clusters, irrespective of redshift and a. 
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Figure 2. Histogram of mean velocity dispersion of the clus- 
ters: a) the total sample; b) "relaxed" clusters. 

5.1. Comparison with other estimates 

20 of the clusters in our sample are also in the ENACS 
survey. Figure 3a shows the comparison between our esti- 
mates and those obtained by Fadda et al. 1996. We found 
a mean difference < (JFadda — osars >— —89 ± 132 that 
indicate that our values of a are on average slightly higher 
than those in Fadda et al. 1996. If we restrict our sample 
to those clusters with at least 30 confirmed members (the 
same restriction is applied for Fadda et al. 1996) we find 
< a 'Fadda - o-sars >= -40 ± 108 which suggest a smaller 
shift and spread (see figure 3b). In both cases we have 
made the comparison assuming the same cluster radius as 
Fadda et al. 1996 (typically smaller than our maximum 
cluster radii). 




1000 1200 



Figure 3. a) Comparison with Fadda et al. 1996 results. The 
solid line corresponds to equal a. b) Same as a) for clusters 
with at least 30 confirmed members in both sample. The solid 
line corresponds to equal a. 



5.2. Cluster mean velocity dispersion vs. richness counts 

Taking into account the methods previously described, 
we have computed the mean velocity dispersions for our 
sample of clusters. We have performed a comparison be- 
tween a and the richness number counts Af as defined in 
the ACO catalog. Since only a small fraction of our cluster 
sample are known X-ray emitters, we have not attempted 
to analyze correlations between our dynamical estimates 
and the X-ray information. 

Figure 4a show the correlation between a and rich- 
ness counts Af taken from the ACO catalogue where no 
clear correlation can be appreciated. A similar result was 
found by Mazure et al. (1996). These authors suggest 
that the very broad relation between Af and a must be 
largely intrinsic. Nevertheless, when we restrict to the 
"relaxed" clusters and "isothermal" distributions (gaus- 
sian velocity distribution and VDP flat or slowly decay- 
ing, see the next section for details) and exclude those 
clusters more strongly affected by projection effects, the 
data suggest some correlation between richness counts and 
cr in the sense that richest clusters tend to have higher 
g. This correlation can be see in Figure 4b, where a lin- 
ear fit has been applied deriving the following relation: 
a = (6.2 ± 2.8)Af + (158 ± 202). 
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Figure 4. a) Correlation between the mean velocity dispersion 
estimate and the richness counts estimated by ACO. b) Same 
as a) for "relaxed" clusters. 



5.3. Velocity dispersion profiles 

The large projected area around clusters in the SARS 
survey allows us to analyze of the dynamics of galaxies in 
the extended halos of clusters. A useful statistical measure 
of the dynamics is the velocity dispersion profile, the veloc- 
ity dispersion at a given radius evaluated by using all the 
galaxies within that radius. The VDP was computed for 
the 29 clusters in our sample with more than 20 confirmed 
members. We have used a step size of 0.5 Mpc hr 1 while 
most VDP are computed up to 4 Mpc h^ 1 and in some 
cases beyond this radius. Many clusters show an irregu- 
lar trend in the inner part (cluster radius < 1 Mpc h^ 1 ), 
this effect could be partially related to the low number of 
galaxies in the inner part of the cluster and may also de- 
pend on the choice of the cluster center which in our case 
correspond to the values provided by Abell. Nevertheless, 
the most interesting aspect of the VDPs is the behavior 
at large distances from the cluster center. For these 29 
clusters we find 19 (14 are "relaxed" clusters) with a flat 
VDP, 5 that present a slowly decaying profile and 5 with 
a rising profile. It should be noted that only 1 of the VDP 
rising clusters was classified as "relaxed" cluster. These re- 
sults are shown in Figure 5 for the 29 objects with reliable 
estimates of VDP. 



7 



E 

Q 
> 



1000 
500 


1000 
500 


1000 
500 


1000 
500 


1000 
500 


1000 
500 



o 



1 1 1 i 1 1 1 

A0080 



II 

A2854 _ 



jiiiilli 



+++ 



A2933 



3 



ii 

A3141 



I I I I I I I I I 



A3194 



i | i i i | i i i | 

tt A3 



i i i I i i i 



1 I 1 1 1 I 
A1750 



I I I I I I I 



A2871a_ 



Iiii 



an 



1 1 1 1 1 



A3111 



I I I I I I I 



A3223a_ 



urn* 



Ml III I I I I 



A3915 



ii 



A2734 



A2871b_ 



I I I I I I I I I 



A3112 



l l l | l l l | l l l 

A3151 



I I I | I I I | I I I 
A3264 



1 1 1 1 1 1 1 1 1 1 1 1 

A3921 



i i i I i i i 



A2799 



I I I | I I I 
A2911 



1 1 1 1 1 1 1 



A3122 



II 

A3153 




I I I I I I I 



A3764 



1 1 1 1 1 1 1 1 1 

A4010 _ 



1 1 1 1 1 1 1 



A2800 J 



II 



A2923 

I i i i I i i i 



A3135 



II 

A3189 1 



ii 

A3809 



I I I I I I I I : I I 



I I I I I I I I 



2 4 60 2 4 60 2 4 60 2 4 60 2 4 6 



Cluster Radii (Mpc/h) 



Figure 5. Velocity dispersion profiles for 29 Abell clusters. 
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In order to allow for a physical comparison between clus- 
ters with different mean velocity dispersions, we have nor- 
malized cluster's radii using r 2 oo (the radius where the 
mean interior cluster overdensity is 200). Assuming a sin- 
gular isothermal profile Carlberg et al. 1997 derive the 
following correlation between r 2 oo and the cluster mean 
velocity dispersion: r 20 o = \qho(z) • ^ e nave followed a 
similar analysis than that proposed by den Hartog and 
Katgert (1996) and Jing and Borncr (1996) consisting in 
the computation of the ratios of a at different distances 
from the cluster center. Both den Hartog and Katgert 
(1996) and Jing and Borner (1996) use the radius in Mpc 
(up to 3 Mpchr x and 1 Mpchr x respectively). We pro- 
pose the use of a normalized radius and four bins for the 
computation of the velocity dispersion estimates a\, er 2 , 
<7 3 and ct 4 : r/r 200 < 1, r/r 2 oo < 2, r/r 200 < 3, r/r 200 < 7 
respectively. The shape of the VDP at different radii can 
be quantified by the ratios <Ji/(?j (=1 for a flat profile). In 
figure 6 we show the distribution of the following ratios: 
ci/cr 2 , 02/(73, 01/03, 01/04 (t 2 /<T4 and the derived mean 
values are 0.93, 1.00, 0.95, 0.96 and 1.00 respectively. It 
can be appreciated that those ratios involving <j\ suggest 
that in the inner bin (up to r/r 20 o = 1) the velocity disper- 
sion is approximately 10% lower than at larger distances. 
This fact can also be appreciated in figure 7 where the to- 
tal sample of clusters has been averaged being each VDP 
profile normalized with the corresponding mean cluster ve- 
locity dispersion. Figure 7 also clearly shows that the av- 
erage VDP for r/r 20 o > 1.5 is nearly flat well beyond the 
cluster virial radius. 
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Figure 6. Velocity dispersion ratios between as computed 
for different normalized bins: a±, 02, 03 and 04 correspond to 
r/r 2 oo < 1, r/r 2 oo < 2, r/r 2 oo < 3, r/r 20 o < 7 respectively. 

6. DISCUSSION AND CONCLUSIONS 

A significant degree of substructure in clusters of galax- 
ies is expected in the hierarchical scenario of structure for- 
mation. This is because of the large time-scale for the 
remnants of the accretion of groups onto clusters in the 
recent past to be erased by dynamical relaxation. On the 
other hand, the identification of clusters in two dimen- 
sions may be strongly biased by spurious systems due to 
projection effects as shown in numerical simulations (van 
Haarlem, Frenk & White 1997). These two issues heavily 
complicate a detailed analysis of the dynamical properties 
of clusters of galaxies. 

The evidence for substructure in rich clusters have been 
extensively explored in different studies of the galaxy dis- 
tribution in cluster fields where the most accepted view 
is the relevant presence of substructure. The X-ray ob- 
servations also contribute to our knowledge of the spatial 
properties of the intra-cluster gas. However, recent analy- 
ses provides conflicting results regarding the isothcrmality 
of the gas or the existence of steeply-declining tempera- 
ture profiles. This issue, and the fact that the degree of 
substructure increases at large distances from cluster cen- 
ters motivated the present study of the outer regions of 
clusters of galaxies. Our work is mainly centered in the 
analysis of the radial velocity distribution of galaxies in 
extended regions of Abell clusters, focusing on the exis- 
tence of gradients in the velocity dispersion profiles. 

We have carried out an analysis of the velocity field of 
galaxies in extended regions up to 7h~ x Mpc from the 
cluster centers. We have applied several methods to re- 
move contamination by projection effects and analyzed the 
presence of sub-clustering. We have obtained suitable es- 
timates of the mean velocity dispersions and its radial de- 
pendence using the ROSTAT routines. Our analysis can 
be compared to Fadda et al. 1996 for a fraction of com- 
mon objects. It is clear from our analysis that the larger 
differences arise in those clusters with more contamina- 
tion and a smaller number of measured redshifts. We also 
find that the correlation between mean velocity dispersion 
a and richness number counts M is strongly affected by 
projection effects. There is some evidence of correlation 
between a and M for a sub-sample restricted to systems 
with no significant degree of contamination. 

From our original sample of 41 Abell clusters we found 
that 40 are real clusters although 4 of these appear as dou- 
ble systems. These results are similar to those found by 
Mazure et al. 1996. These authors found that almost all 
ACO clusters with richness class 1 or greater correspond 
to real systems in the redshift space and about 10% of the 
ACO clusters appears to be the result of a superposition of 
two similar poorer systems. Beside the double systems we 
also found that 7 of the clusters in our sample are subject 
to serious projection effects. The fraction of clusters with 
a high degree of contamination in our sample compares 
well with the results of such effects in the mock catalogs 
from the numerical simulations of van Haarlem, Frenk & 
White 1997) where 1/3 of Abell-type clusters are expected 
to arise from projection of groups along the line of sight. 
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From the resulting sample of 44 clusters, four are poorly 
denned and more data are needed in order to better estab- 
lish their properties. In spite of projection effects, 28 of 
the 44 clusters are well defined in redshift space and have 
a velocity distribution consistent with a gaussian function. 

Our results show that the average VDP is flat at large 
distances from the cluster center. This behavior is found 
for 19 clusters (65% of the 29 with VDP estimates) indi- 
cating that an isothermal hypothesis can be assumed even 
at radii well beyond the virial radius. Nevertheless, we 
found that on average, the normalized velocity dispersion 
is about 10% smaller in the inner region of the clusters 
(r/r2oa < 1). A possible interpretation for the decay of 
VDPs in central regions can be related to the morpholog- 
ical segregation in clusters. Fadda et al. 1996 found kine- 
matical segregation in the sence that early type galaxies 
show smaller values of a than late types and Mazure et 
al. 1996 found that the brightest cluster galaxies (typi- 
cally of early type morphology) move slower than other 
galaxies and Ramirez et al. 1998 found that the differ- 
ences between the velocity distributions of elliptical and 
spiral galaxies are associated with the shape of their orbit 
families. Since early type objects are dominat within r2oo 
the results shown in figure 7 are to be expected. 



Figure 7. Mean VDP for the total sample of clusters with 
VDP estimate. Each cluster has been normalized using the 
corresponding mean a. 

The shape of VDP profiles are of fundamental impor- 
tance for their implications on cluster properties and cos- 
mology since detailed theoretical predictions from differ- 
ent cosmological scenarios could be used to set restric- 
tions to current models of structure formation. Jing and 
Bonier 1996 analysis of VDPs of clusters for several cos- 
mological models show an average decline of VDPs with 
the distance to cluster centers. Nevertheless, this analy- 
sis was restricted to the very inner region of the clusters 
(< lh~ 1 Mpc. Therefore, new numerical simulations must 
be analyzed to test our findings of flat VDPs at very large 
distances from the cluster center. 

Under the assumption that clusters are in global dynam- 
ical equilibrium (even beyond the virialization radius) our 
results can be compared to temperature radial distribu- 
tions derived from the X-ray emission of the intra-cluster 
medium. Flat VDP at large clustercentric distances may 
shed light on the recent controversy on the nature, either 
flat or declining, of intracluster temperature radial profiles 
(see White 2000). 
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